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S3 ■ Abstract 

We study the spectral distribution of the threshold network model. The results 
contain an explicit description and its asymptotic behaviour. 
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^ 1 Introduction 



The threshold network model ^„(X, 6), where X is a random variable, n > 2 is an integer 
O I and ^ G M is a constant called a threshold, is a random graph on the vertex set V = 
{1,2,..., n} obtained as follows: let Xi, X2, . . . , X„ be independent copies of X and draw 
an edge between two distinct vertices i,jEV if Xi + Xj > 6. In other words, Qn{X, 6) is 
. specified by the random adjacency matrix A = {Aij) defined by 

rS 

1 0, otherwise, 

where Is denotes the indicator function of a set I? . As a small variant one may allow 
self-loops, see e.g., |4j. In this case the threshold network model is denoted by 6*), 
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where two vertices i,j & V (possibly i = j) are connected if Xj + Xj > 9. The adjacency 
matrix A = (Aij) is given by 

Aij = I^g^^){Xi + Xj), i,jeV. 

The threshold network model has been extensively studied as a reasonable candidate 
model of real world complex graphs (networks), which are often characterized by small 
diameters, high clustering, and power-law (scale- free) degree distributions (HEIEQ]. In 
fact, the threshold network model belongs to the so-called hidden variable models [51122] 
and is known for being capable of generating scale-free networks. Their mean behavior 
[3l|5l|8l|9l[T5l[2Tl|22] and limit theorems [Tt fTTHTS] for the degree, the clustering coefficients, 
the number of subgraphs, and the average distance have been analyzed. See also [^[TTl^ 
[lllIinillTj for related works. 

Spectral properties of the threshold network model are also of interest. As a simple 
case, the binary threshold model appears in [23]. The strong law of large numbers and 
central limit theorem for the rank of the adjacency matrix of the model with self- loops 
are given by [1] . Eigenvalues and eigenvectors of the Laplacian matrix of the model have 
been studied [ISIIII]. For general results of spectral analysis of graphs see e.g. Hora- 
Obata [To]. The main purpose of this paper is to study the spectral distribution of 
the threshold network model. Our result covers the preceding study of the rank of the 
adjacency matrix. 

This paper is organized as follows: In Section 2 we recall the hierarchical structure of 
the threshold network model and derive the spectral distribution of each sample graph 
(threshold graph). In Section 3 we obtain similar results for the threshold network model 
which admits self-loops. In Section 4 we derive some asymptotic behaviors of the spectral 
distributions and in Section 5 we give a simple example called the binary threshold model. 

2 Spectra of threshold graphs 

Each sample graph G G Qn{X,6) has a hierarchical structure described by the so-called 
creation sequence, introduced by Hagberg-Schult-Swart [9]. Here we adopt a variant by 
Diaconis-Holmes-Janson [6] . Each G being determined by the values of random variables 
Xi, X2, . . . , Xn, we arrange them in increasing order: X(i) < X(2) < ■■■ < X(„). If 
X(i) -|- > 9, we have 

6 < X(i) + X(„) < X(2) + < ■ ■ ■ < X(„_i) + X(„), 

which means that the vertex corresponding to X(„) is connected with the n — 1 other 
vertices. Otherwise, we have 

& > -^(1) + > ■ ■ • > X(i) + X(3) > X(i) + X(2), 

which means that the vertex corresponding to X(i) is isolatedD We set = 1 or s„ = 
according as the former case or the latter occurs. Then, according to the case we remove 
the random variable X(„) or X(i), we continue similar procedure to define s„_i, . . . ,S2. 
Finally, we set Si = S2 and obtain a {0, l}-sequence {si, S2, . . . , which is called the 
creation sequence of G and is denoted by Sq- 
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Given a creation sequence Sq let ki and k denote the number of consecutive bits of 
I's and O's, respectively, as follows: 



ki h k'z h km Im 

Sg = {C^, oT^, C^, oT^To, . . . , C^, C^}- (1) 

It may happen that A;i = or = 0, but we have k2, ■ ■ ■ , km, h, . . . , Im-i > 1 and m > 1. 
Moreover, by definition we have two cases: (a) ki = (equivalently Si = 0) and li > 2; 
(b) ki >2 (equivalently, si = 1). 

For example, if Sq = {1, 1,0,0, 1,0, 1,0} then ki = 2, li = 2, k2 = 1, h = 1, ^3 = 
1, /a = 1 and Fig. [T] shows the shape of G. 
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Figure 1: A threshold graph G corresponding to Sg = {1,1,0,0, 1,0, 1,0} 



The creation sequence Sg gives rise to a partition of the vertex set: 

m m 



i=l 1=1 



The subgraph induced by v}^^ is the complete graph of ki vertices, and that induced 
by V^*"'^'' is the null graph of k vertices. Moreover, every vertex in v}^^ (resp. V^^^) is 
connected (resp. disconnected) with all vertices in 

In general, a graph possessing the above hierarchical structure is called a threshold graph 



Theorem 1. Let G be a thereshold graph with a creation sequence Sg = {si = S2, S3, ... , s„}. 
Define ki and k as in ([T]) and set 



m 



Gn{-1) = J2k^-{nl-l)- /{l}(Si), Gn{0) = J] - (m - 1). (2) 



i=l i=l 

Then the spectral distribution of G is given by 



/in(G) = ^^^^5_i + ^5o + -V5A,, J = 2(m-l) + /{i|(.0, (3) 

n mm < ^ ■' ^ ' 



where {Xj} exhausts the eigenvalues of the matrix: 
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(5) 



for Si = (equivalently, ki = 0). Moreover, any Xj in ([3]) differs from and —1. 

Proof. Let lij denote the ^ x j matrix consisting of only 1, Ojj the i x j zero matrix, 
li the i X i identity matrix, and Ij^j = Ij^j — /j. By the hierarchical structure mentioned 
above, the adjacency matrix of G is represented in the form: 



A 



G 



0; I 
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)^m — 1 


5^m — 1 


^^m?^m — 2 






Ofc I 
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The adjacency matrix A acts on C" from the left. We define subspaces of C" by 



V^.(O) 

K^(0) 




, 1 < i < m, 



6 + 6 + ■ ■ ■ + a. = 

Vi + V2-\ \-Vh = ^ 1 < i < m - I, 
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where 
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Since Ag- acts on Vi{—1) as the scalar operator with —1, it possesses the eigenvalues — 1 
with multiplicity at least 



m 



5^dimy,(-l) = 5^(A:,-l) = 5^A:,- 

i=l i=l i=l 

ii ki>2 (i.e., Si = 1), and 

m mm 

dim Vi{-1) = J2iki - 1) = X A;, - (m - 1) 



i=2 



i=2 



i=2 



if /ci = (i.e., si = 0). In any case, the multiphcity is at least C„(— 1) defined in (|2]). 
Similarly, acting on 1^(0) as a scalar operator with 0, Aq possesses the eigenvalues with 
multiplicity at least C„(0). 

Let W be the orthogonal complement to (V^(— 1) © ^.(0)). The matrix represen- 
tation of Ag on W with respect to the basis 



0, 







1 < i < m, and Wi 



1 < i < m — 1. 



is given by (|4]) or by ([5]) according as fci > 2 or fci = 0. Then, one may verify easily the 
eigenvalues of the matrices (jl]) and are different from —1 nor 0. □ 

Remark After simple calculation we see that the eigenvalues Ai, . . . , Aj in ([3]) are ob- 
tained from the characteristic equations 



M(A) = 0, 



where 



M(A) = det 
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if fci > 2 (i.e., si = 1), and 



M(A) = det 



km 1 A 
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if fci = (i.e., Si = 0). Simple calculation shows that 
M(-l) = 

and 



ki - ■ -km-h - ■ ■ Im-i, if fci > 2 (i.e., si = 1) 

k2 - ■ ■ km - {k - I) ■ h- ■ ■ Im-i, otherwise. 



M(0) 



{ki - 1) ■ k2- ■ ■ km ■ h- ■ ■ Im-i, if fci > 2 (i.e., si = 1), 
k2 ■ ■ ■ km ■ h ■ ■ ■ lm-ii otherwise, 

from which we see also that do not contain —1 or 0. 



3 Spectra of threshold graphs with self-loops 

The idea of a creation sequence in Section [2] can be applied to the threshold network 
model which allows self-loops. With each G G Qn{X,9) we associate a creation sequence 
Sg = §25 • • • 5 Sn} as follows: if X(i) + X(„) > 6, we have 

which implies that the vertex corresponding to X(„) is connected with the n — 1 other 
vertices and has a self-loop. Otherwise, 

d > -'^(l) + X(^ri) > ■ ■ ■ > -^(1) + -^(3) > -^(1) + -^(2) > -^(1) + -^(1) , 

which means that the vertex corresponding to X(i) is isolated and has no self-loopsD We 
set s„ = 1 or s„ = according as the former case or the latter occurs. Then, according 
to the case we remove the random variable Xi^n) or we continue similar procedure to 
define Sn-i, . . . ,S2- Finally, letting X(^,) be the last remained random variable, set si = 1 
if X* > 6/2 and Si = otherwise. In this case G is called a threshold graph with self-loops 
associated with a creation sequence S = {si, S2, . . . , We note that if Sj = 1 the 
corresponding vertex has a self-loop and otherwise no self-loop. 

Given a creation sequence S = {si, §2, ■ ■ ■ , Sn} we define kj and Ij as in ([T]). It may 
happen that ki = and Im = 0, but k2, . . . , km, h, ■ ■ ■ , Im-i > 1 and m > 1. The adjacency 
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matrix of G is of the form: 
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(6) 

Repeating a similar argument as in Theorem [U we come to the following 

Theorem 2. Lei G be a threshold graph with self-loops associated with a creation sequence 
S = {si,S2, . . . ,Sn} and its adjacency matrix given as in Set 

C„(0) = n-2(m-l)-/|i}(5i). (7) 

Then the spectral distribution of G is given by 

MG) = ^ 5o + - 5^ 5a, , J = 2{m - 1) + (8) 

where {Xj} exhaust the eigenvalues of 
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for si = (i.e., ki = 0). Moreover, any Xj in ([8]) differs from 0. 
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Remark The eigenvalues Ai, 
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for si = (i.e., ki = 0). 



4 Limit theorems 

In this section we discuss asymptotic behaviors of the spectral distributions obtained in 
the previous sections. 

We first consider the case where the distribution of X is discrete and given by 

oo 

r{x = i)=p,, 1 = 0,1,..., X]P' = i- 

i=0 

Let m > 1 be a fixed integer. Take a particlar threshold 6 = 2m — 1 and assume that 
Pi > for i = 0,1, . . . , 2m — 1. It follows from the strong law of large numbers that 

k = : Xj = m - i}, i = l,...,m, 
ki = 'ilij : Xj = m - 1 + i}, i = 1, . . . ,m - 1, 
kra = tl{j : Xj > 2m - 1}, 
h = ki = 0, i > m + 1, 
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for large n almost surely. Moreover, denoting by F the distribution function of X, we 
have 

lim — > /j = F{m — 1) a.s. lim — N fcj = 1 — F(m — 1) a.s. 

n— !>oo n — ^ n— )-oo n — ' 

With these observation we easily obtain the following 

Theorem 3. Notations and assumptions being as above, the spectral distributions of 
Qn{X^ 2m — 1) verifies 

lim = (1 - F{m - 1)) ■ 5_i + F{m - I) ■ 5o a.s. 

n—^oo 

Similarly, the spectral distributions o/^„(X, 2m — 1) verifies 

lim Jln{G) = 6o a.s. 

n—^oo 

Remark Similar results hold when the distribution of X is discrete and F has only finite 
number of jumps in (— oo,^/2] or [6/2, oo). But no simple description is known for a 
general case. 

Next we consider the case where the distribution of X is continuous. As is stated by 
Bose-Sen |3] implicitly, if the distribution of X is continuous and symmetric around 0, 
then the distribution of zero and one entries in the creation sequence S of each graph 
generated by ^„(X, 0) is the same as the distribution of a sequence of i.i.d. Bernoulli 
random variables {l^j}j=i,2,...,n with success probability 1/2, that is, 

P(y. = 0) = F(Yi = 1) = 1/2, for z = 1,2,.. . ,n. 

This means that S = {si, S2, . . . , s„} = {Yi,Y2, . . . , 1^}. Recall that the creation sequence 
5* of each graph generated by Qn{^, 0) is always satisfied with si = S2. Then we observe 

that S = {si = S2, S3, ... , Sn} = {Y2, ¥2, 13, • • • , K}, where {Fj}i=2,3,...,n be the sequence 
of i.i.d. Bernoulli random variables with success probability 1/2, similarly. 

Taking the above consideration into account, we obtain asymptotic behaviors of coef- 
ficients of point measures on —1 and appearing in fin{G) and Jln{G). 

Theorem 4. Assume that the distribution of X is continuous and symmetric around 0. 
Define C„(— 1), C„(0) and C„(0) as in ([2]) and ((Tj). Then we have 

(1) lim Cni-l)/n = lim C„(0)/n = 1/4 a.s. 

ra— >oo ra— >oo 

(2) v^(C„(-l)/n- 1/4) A^(0,l/4) and (C„(0)/n - 1/4) ^ A^(0, 1/4) as n ^ 
00. 

(3) lim CniO)/n = 1/2 a.s. 

n—^oo 

(4) ^{CniO)/n - 1/2) iV(0, 1/4) as n 00. 
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Proof. Note the following relations: 
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i=l i=l 



We then easily check that 

1 n 



E[a(-1)]=E[C„,_,, 2 4 
and 

E[C„(0)] = ^ . 

Applying a similar argument as in [H Theorem 1], we have the assertion. □ 

When the distribution of X is continuous and symmetric around 9/2, we can obtain 
similar results for ^n(X, 6) and Qn{,X, 6) by straightforward modification. Study covering 
a more general situation is now in progress. 



5 Binary threshold model 

In this section we give a simple example. The threshold network model defined by 
Bernoulh trials Xi,X2, . . . ,X„ with success probability p, i.e., < = 1) = p < 1, 

and a threshold < 6* < 1 is called the binary threshold model and is denoted by Qnip)- 
For G G Qnip) the partition of the vertex set V is given by 

\/ = \/Wuy(°), = {z; = 1}, = {z; =0}. 

Theorem 5. For G E Qn{p) we set \V^^^\ = k and \ V^^^\ = I. Then the spectral distribu- 
tion of G is given by 

A;-1 /-I 1 , 1 

IJ-k,! = ^-1 H Oo + - dx+ + - Ox_ , 

n n n n 
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where 

k-i±^ik-iy + m 

^± = 2 ■ 

Proof. We need only to apply Theorem [1] with li = I, I2 = ki = 0, k2 = k and m = 2. In 
this case, (EI) becomes 

\k-l I 
k 

of which the eigenvalues are A-t in (|9]). □ 
Corollary 1. Let Hn{G) be the spectral distribution of G E Gnip)- Then we have 

lim /i„(G') = p ■ 6-1 + {1 - p) ■ 60 a.s. 

n— >oo 

Proof. By the strong law of large numbers, see also Theorem [3 □ 

As for the the mean spectral distribution we have 

Theorem 6. The mean spectral distribution of the binary threshold model Qn{p) is given 
by 

fin= (p - H + fl - P- H 



k=0 

where 



I — n \ / 



(10) 



Proof. Since 



k-l±^{k-iy + Ak{n-k) 
X±{k) = , A; = 0,1, 



P(|\/(i)| = k, |\/(°)| = /) = - pY, k + l = n, 



the mean spectral distribution is given by 



k=0 ^ ^ 



Then ( ITOl) follows from Theorem [5] by direct computation. □ 

Corollary 2. Let fin be mean spectral distribution of the binary threshold model Qn{p)- 
Then we have 

lim /i„ = p ■ 5_i + (1 - p) ■ 5o ■ 

n— )-oo 
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